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ABSTRACT 

In weakly-collisional plasmas such as the intracluster medium (ICM), heat and momentum 
transport become anisotropic with respect to the local magnetic field direction. Anisotropic 
heat conduction causes the slow magnetosonic wave to become buoyantly unstable to the 
magnetothermal instability (MTI) when the temperature increases in the direction of gravity 
and to the heat-flux-driven buoyancy instability (HBI) when the temperature decreases in the 
direction of gravity. The local changes in magnetic field strength that attend these instabilities 
cause pressure anisotropics that viscously damp motions parallel to the magnetic field. In this 
paper we employ a linear stability analysis to elucidate the effects of anisotropic viscosity (i.e. 
Braginskii pressure anisotropy) on the MTI and HBI. By stifling the convergence/divergence 
of magnetic field lines, pressure anisotropy significantly affects how the ICM interacts with 
the temperature gradient. Instabilities which depend upon the convergence/divergence of mag- 
netic field lines to generate unstable buoyant motions (the HBI) are suppressed over much 
of the wavenumber space, whereas those which are otherwise impeded by field-line conver- 
gence/divergence (the MTI) are strengthened. As a result, the wavenumbers at which the HBI 
survives largely unsuppressed in the ICM have parallel components too small to rigorously be 
considered local. This is particularly true as the magnetic field becomes more and more or- 
thogonal to the temperature gradient. The field-line insulation found by recent numerical sim- 
ulations to be a nonlinear consequence of the standard HBI might therefore be attenuated. In 
contrast, the fastest-growing MTI modes are unaffected by anisotropic viscosity. However, we 
find that anisotropic viscosity couples slow and Alfven waves in such a way as to buoyantly 
destabilise Alfvenic fluctuations when the temperature increases in the direction of gravity. 
Consequently, many wavenumbers previously considered MTI-stable or slow-growing are in 
fact maximally unstable. We discuss the physical interpretation of these instabilities in detail. 

Key words: conduction - instabilities - magnetic fields - MHD - plasmas - galaxies: clus- 
ters: intracluster medium. 



1 INTRODUCTION 

The intracluster medium (ICM) is stratified, not only in pressure, 
but also in entropy. Until recently, the latter was thought to be of 
paramount importance to the ICM's dynamical stability. The rea- 
son for this is easy to understand. In an atmosphere where entropy 
increases upwards, an upward adiabatic displacement of a fluid el- 
ement leaves the element cooler than its surroundings. A cool el- 
ement is denser (because of local pressure balance), so the buoy- 
ancy force is restoring. If, on the other hand, the entropy were to 
decrease upwards, an upward adiabatic displacement produces a 
fluid element that is warmer than its surroundings and so there is 
no restoring buoyancy force, t he fluid element conti nues to rise, and 
convective instability ensues dS chwarzschildll 1 95 81) . Careful obser- 
vations have shown that the ICM generically has a positive entropy 
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gradient (e.g. ICavagnolo et alj|2009l) and so, by this reasoning, is 
convectively stable. 

In fact, matters are not this simple. The ICM is not a con- 
ventional fluid. Firstly, it is magnetized. Secondly, particle-particle 
collisions are rare. The conductive flow of heat consequently be- 
comes strongly anisotropic with respect to the local magnetic field 
direction, since collisional energy exchange by (predominantly) the 
electrons can occur much more readily along the magnetic field 
than across it. The heat is then restricted to being channeled along 
magnetic lines of force. When the conduction timescale is shorter 
than any other dynamical timescale in the system, magnetic field 
lines become isotherms. 

Such a radical change in the thermal behaviour of a weakly 
collisional plasma profoundly alters its stability properties. Draw- 
ing on the powerfu l analogy betw een angular momentum and en- 
tropy stratification. iBalbusl 112000) argued that the temperature gra- 
dient takes precedence over the entropy gradient in determining 
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the convective stability of a weakly collisional plasma. He sup- 
ported this conjecture with a linear stability analysis that proved 
the existence of th e magnetothermal instability (hereafter, MTI; 
Balbus 200(31 l200lh . which is triggered in regions where the tem- 
perature (as opposed to the entropy) inc reases in the di r ection 
of gra vity. Subsequent numerica l work by IParrish & Stonel (I2005L 
2007) and McCourt et al.l l l2010h demonstrated the efficacy of the 
MTI and extended it into the non-linear regime. Numerical stud- 
ies of the effect of the MTI on the evolution of the outer re- 
gions of non-isothermal galaxy clusters, where the temperature de- 
crease s with distance from the central core, followed soon there- 
after ( IParrish. Stone & Lemas ter 2008). 

In contrast, the inner ~200 kpc or so of non-isothermal 
clusters are characterised by outwardly increasing temperature 
profile s (e.g. IPiffaretti et alj|2005l: IVikhlinin et al.ll2005h . lOuataertl 
(2008) showed that the ICM in these inverted temperature profiles 
is also buoyantly unstable to a heat-flux-driven instability, now 
referred to as the HBI. The HBI arises because perturbed fluid 
elements are heated/cooled by a background heat flux in such a 
way as to become buoyantly unstable. Recently there has been a 
surge of numerical efforts to understand the non-linear evolution 
of the HBI and its implications for the so-c alled 'cooling-flow 
probl e m' exhibited by cool-core clusters (jP arrish & Ouataert 
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In this pa per, we extend the work of iBalbusi ( 120001) and 
lOuataerj j2008l) to include the effects of pressure anisotropy (i.e. 
anisotropic viscosity). We are motivated by the simple consider- 
ation that one cannot self-consistently take the limit of fast ther- 
mal conduction along magnetic field lines while simultaneously 
neglecting differences between the thermal pressure parallel and 
perpendicular to the magnetic field direction. While the dynami- 
cal effects of pressure anisotropy occur on a timescale longer than 
conduction, in weakly collisional plasmas such as the ICM this 
timescale is nevertheless still shorter than (or at least as short as) 
the dynamical timescale and, as a consequence, the MTI and HBI 
growth times. Our principal result is that, by stifling the conver- 
gence/divergence of magnetic field lines, pressure anisotropy sig- 
nificantly affects how the plasma in the ICM interacts with the 
temperature gradient. Instabilities which depend upon the conver- 
gence/divergence of magnetic field lines to generate unstable buoy- 
ant motions (the HBI) are suppressed over much of the wavenum- 
ber space, whereas those which are otherwise impeded by field-line 
convergence/divergence (the MTI) are strengthened. We also com- 
ment on how pressure anisotropy affec ts the heat-flux-driven bu oy- 
ancy overstability recently found bv lBalbus & Reynolds! d2010h . 

The paper is organised as follows. In Section [2] we motivate 
the inclusion of pressure anisotropy in our picture of weakly colli- 
sional buoyancy instabilities and provide a qualitative discussion 
of its effects. Readers not interested in the mathematical details 
may read this Section and proceed immediately to the conclusions 
(§2). In Section[3] we formulate the problem by presenting the ba- 
sic equations, linearising them, and deriving the dispersion relation 
governing small perturbations about a simple equilibrium state. The 
solutions of this dispersion relation are examined in Section|4] We 
conclude in Section [5] with a summary of our results and a brief 
discussion of their implications for the structure and evolution of 
the ICM. 



2 PHYSICAL MOTIVATION AND THEORETICAL 
EXPECTATIONS 

The HBI and MTI operate most efficiently at sufficiently small 
wavelengths for which the conduction rate is much greater than 
the local dynamical frequency, i.e. cj con d S> £Jd yn = io/H) 1 ^ 2 = 
Vth/H, where g is the gravitational acceleration, H is the thermal- 
pressure scale-height of the plasma, and v t h is the thermal speed 
of the ions. This precludes the usual buoyant restoring force, which 
would otherwise result in Brunt- Vaisala oscillations, by ensuring 
magnetically-tethered fluid elements communicate fhermodynami- 
cally much faster with one another than they do with the ambient 
medium. In weakly collisional environments such as the ICM, this 
timescale separation is satisfied for a wide range of wavenumbers 
satisfying fc|| (A m f p H) 1//2 S> (m /mi) 1//4 ~ 0.1, where fc| is the 
wavenumber along the magnetic field, A m f p = Uth/^i is the par- 
ticle mean free path between collisions, and v\ is the ion-ion col- 
lision frequency. Typical values of i//A m f p in the ICM decrease 
outwards from ~10 3 to ~10 2 in the cool-cores of non-isothermal 
clusters, and from ~10 2 to ~10 beyond the cooling radius out to 
~1 Mpc. It is important to note that, while the HBI and MTI owe 
their existence to rapid conduction along magnetic field lines, the 
unstable perturbations themselves only grow at a rate ~ojd yn . This 
is because the free energy required to drive the instabilities is ex- 
tracted from the background temperature gradient, which is set by 
macroscale processes, at a rate determined by gravity. 

There is however another timescale that ought to be consid- 
ered. In a magnetized plasma, any change in magnetic field strength 
must be accompanied by a corresponding change in the perpendic- 
ular gas pressure, since the first adiabatic invariant for each particle 
is conserved on timescales much longer than the inverse of the ion 
cy clotron frequency (which is extremely short in the ICM; see §2.2 
of lKunz et alj|201 3 ). The resulting p ressure anisotropy is the phys- 
ical effect behind what is known as iBragin skii ( 1965) viscosity - 
the restriction of the viscous damping (to dominant order in the Lar- 
mor radius expansion) to the motions and gradients parallel to the 
magnetic field[jThis in turn implies that perturbations in magnetic 
field strength are erased at the viscous damping rate oj v i sc (see eq. 
I23t . While cj v i sc is smaller than cj con d by a factor ~10 (see end of 
Section |33l l, and so it may be tempting to ignore viscosity relative 
to conduction, oj v i sc is in fact much greater than the growth rates of 
the MTI and HBI at most wavelengths at which the latter are usu- 
ally thought to operate in the ICM. Here we provide a qualitative 
discussion of these effects. 

The HBI relies on the presence of a background heat flux, 
which may be tapped into by the convergence and divergence of 
conducting magnetic field lines. Downwardly-displaced fluid ele- 
ments find themselves in regions where field lines diverge; they 
are conductively cooled via the background heat flux, lose energy, 
and sink further down in the gravitational potential. As they sink, 
the local field lines diverge further and an instability ensues. By 
contrast, an upwardly-displaced fluid element gains energy from 
the converging heat flux and thus buoyantly rises. Braginskii vis- 
cosity hinders the HBI by damping perturbations to the magnetic 
field strength and thereby preventing convergence and divergence 
of field lines (see Fig. [TJ. If oj v i sc ^> ojdyn, the convergence (di- 
vergence) of field lines responsible for the HBI is wiped away 
faster than upwardly (downwardly) displaced fluid elements can 



1 Pressure anisotrop y also leads to microscale plasma instabilities (e.g. see 
Schekochihin et al. 2005, and references therein) but here we will consider 
perturbations around equilibria that do not trigger those. 



Dynamical stability of the ICM 3 



HBI 



buoyancy 
force 



MTI 



with Braginskii viscosity jj Braginskii 



g vr 



z 




*B = B: 



Figure 1. HBI subject to Braginskii viscosity. The plasma is threaded by 
a vertical magnetic field (dashed lines) and has a background heat flux in 
the — z direction. A perturbation (black arrows) with non-zero k x and k z 
modifies the field lines as illustrated (black curves). The heat flux, forced to 
follow the perturbed field lines, converges and diverges, leading to heating 
and cooling of the plasma. For a plasma with dT/dz > 0, a downwardly- 
displaced fluid element loses energy, causing it to sink deeper in the gravi- 
tational field (and vice-versa for an upwardly-displaced fluid element). The 
buoyancy force responsible for this behaviour is denoted by the white solid 
arrows. The pressure anisotropy, which is generated by motions along the 
background field lines, contributes a Braginskii (viscous) force (denoted by 
the grey solid arrows) that impedes this motion. For wavelengths such that 
Wvisc 2> <^dyn> tne two forces become nearly equal and opposite and the 
convergence (divergence) of field lines responsible for the HBI is wiped 
away faster than upwardly (downwardly) displaced fluid elements can take 
advantage of the increased (decreased) heating. 



take advantage of the increased (decreased) heating. In fact, as 
we show in Section |4~T1 the buoyancy and viscous forces become 
nearly equal and opposite when the background field is vertical and 
^visc 3> Wdyn- If the fluid element were to rise buoyantly, it would 
locally increase the magnetic field strength and generate a pressure 
anisotropy, which would cause a viscous stress that damps the ver- 
tical motion and halts the HBI. Pressure anisotropy can be therefore 
thought of as providing an effective tension that 'tethers' a buoyant 
fluid element to its original location, preventing it from rising. The 
only HBI modes to evade strong suppression are those which have 
wavenumbers satisfying o; v isc < i^dyn < w con d; these modes are 
not the same as those usually thought to be the fastest growing. 

Matters are slightly more complicated with the MTI, which 
owes its existence to the alignment of isothermal magnetic field 
lines with the background temperature gradient. It is this align- 
ment that allows a downwardly (upwardly) displaced fluid element 
always to be cooler (warmer) than the surroundings it is passing 
through, even as its temperature rises (falls). As the separation be- 
tween magnetically-connected fluid elements grows, they take the 
magnetic field lines with them, aligning these heat conduits ever 
more parallel to the background temperature gradient and reinforc- 
ing fluid displacements. There is a weak preference for perturba- 
tions whose wavevectors are aligned with the background magnetic 
field; otherwise the consequent convergence (divergence) of any 
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Figure 2. (Top panel) A stable MTI mode, with Braginskii viscosity ig- 
nored, viewed from the side (solid lines) and from above (dots and crosses). 
The background magnetic field B (dashed lines) makes an angle of 45 de- 
grees with respect to both gravity and the temperature gradient. The field 
is perturbed with a wavevector that has equal components in the plane 
perpendicular to gravity, as indicated by the vector k. The corresponding 
eigenvector has SB X = 6B y = 0. This mode is MTI stable because the 
destabilising transfer of entropy from one fluid to another is exactly offset 
by the stabilising exchange of entropy by the convergence/divergence of 
the background heat flux. (Bottom panel) The same mode becomes unsta- 
ble when Braginskii viscosity is self-consistently included. Rapid parallel 
viscous damping effectively orients the perturbed magnetic field nearly per- 
pendicular to the background field, thereby precluding any stabilisation by 
the background heat flux. The perturbed magnetic field has components in 
all three directions in order to satisfy the divergence-free constraint (namely, 
5B Z « 8 By = -5B X ). 



background heat flux would heat (cool) downwardly (upwardly) 
displaced fluid elements (the exact opposite of what happens with 
the HBI), thereby undermining the destabilising upward entropy 
transfer between magnetically-connected fluid elements. 

One effect of pressure anisotropy is to reinforce this prefer- 
ence by suppressing motions along field lines. Another potentially 
more important effect is to destabilise many wave modes that were 
previously thought to be stable to the MTI. One example is shown 
in Fig. [2] The panel on the top exhibits a mode that is stable to 
the standard MTI: the destabilising effect of generating a heat flux 
along field lines is exactly offset by the stabilising effect of con- 
verging/diverging background heat flux. When Braginskii viscos- 
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ity is included and u v isc 3> ^dyn. the same exact mode is unsta- 
ble with a growth rate ~cj c iyn (bottom panel). Any motion along 
the background magnetic field is rapidly damped and so the back- 
ground heat flux cannot interfere with the MTI. In regions where 
the background heat flux converges in the x-z plane it diverges in 
the x-y plane, giving no net heat extraction. In effect, rapid Bragin- 
skii damping effectively endows slow-mode perturbations (which 
are subject to buoyancy forces) with Alfvenic characteristics (i.e. 
perturbed magnetic fields and velocities that are predominantly ori- 
ented perpendicular to the background magnetic field). 

In either case, the fundamental point is that Braginskii viscos- 
ity always suppresses perturbations that acquire free energy from 
the temperature gradient via a background heat flux (see eq. 1341 ). 
This is generally bad for the HBI and good for the MTI. In the next 
Section we formalise these qualitative expectations. 



3 FORMULATION OF THE PROBLEM 
3.1 Basic Equations 

The fundamental equations of motion are the continuity equation 



dp 
At 



-pV-v, 
the momentum equation 



£ = -±v. (p + i£ 

dt p V 87i " 



BB 

4.7T 



+ 9, 



and the magnetic induction equation 
dB 



dt 



B Vv - BV v, 



(1) 



(2) 



(3) 



where p is the mass density, v the velocity, B the magnetic field, 
and g the gravitational acceleration; d/dt = d/dt + ivV is the 
convective (Lagrangian) derivative. 

In the momentum equation (2), the sum of the ion and electron 
pressures 



■p±\- (p± —p\\) bb 



(4) 



is a diagonal tensor whose components perpendicular and par- 
allel (p||) to the background magnetic field direction b = B/B 
are in general different. Differences between the perpendicular and 
parallel pressure in a magnetized plasma arise from the conserva- 
tion of the first and second adiabatic invariants for each particle on 
timescales much greater than the inverse of the cyclotron frequency. 
When the magnetic field strengt h and/or the density change, p± 
and p| | change in different ways {Ch ew, Goldberger & Lowlll956l) . 
For example, conservation of the first adiabatic invariant p — 
mvj_/2B implies that an increase in magnetic field strength must 
be accompanied by a corresponding increase in the perpendicular 
pressure, p±/B ~ const. 

When the collision frequency is larger than the rates of change 
of all fields (i.e. v 3> d/dt) - a condition easily satisfied for buoy- 
ancy instabilities in the ICM - it is straightfor ward to obtain an 
equat ion for the pressure anisotropy (e.g. see ISchekochihin et all 
l20ld for a simple derivation): 

<--«-*5-?w-?( tt ' v -H- ,5 > 



where p — (2/3)pj_ + (l/3)pn is the total plasma pressure^ 
To obtain the final equality, we have used equations (TJ and 10 
to express the rates of change of the magnetic field strength and 
density in terms of velocity gradients. This is referred to as the 
iBraginskiil Jl965t) anisotropic viscosity. The ion contribution to the 
Braginskii viscosity dominates over that of the electrons by a factor 
-(mi/me) 1 / 2 . 

We will also require an internal energy equation, 



3 d p 
2 P dt hl p^ 

where 



= (P± ~P\\ 



d l„ B 



2/3 



V- (bQ 



= - X cb-VT 



(6) 



(7) 



is the collisional heat flux, T is the temperature, and \c is the ther- 
mal conductivity of the electrons 



6 x 10~ 7 T 5/2 ergs cm" 



(8) 



the ther mal conducti vity of the ions is a factor ^(me/mi) 1 ' 2 
smaller JSpitzeJl9 62)Fl Equation expresses the fact that, in the 
presence of a magnetic field, heat is restricted to flow along mag- 
netic lines of force when the particle gyroradi us is much small er 
than the mean free path between collisions (e.g. Braginskii 19651) . 



3.2 Background equilibrium and perturbations 

For simplicity, we consider a plasma stratified in both density and 
temperature in the presence of a uniform gravitational field in the 
vertical direction, g = — gz. The plasma is not self-gravitating, so 
that g is a specified function of position. Without loss of generality, 
the magnetic field is oriented along b — b x x + b z z. We take the 
background pressure to be isotropic and Ti = T c = T, so that 
Pi — p c = p/2. We further assume that the ratio of the ion thermal 
and magnetic pressures is large: 



8-7rpi 



<4 



>i, 



(9) 



where v t h = (p/p) 1/2 = (2k B T/mi) 1/2 and v A = B/(4ttp) 1/2 
are the thermal and Alfven speeds of the ions, respectively. Obser- 
vations of synchrotron radiation, inverse Compton emission, and 
Faraday rotation suggest a plasma f3 parameter that ranges from 
~10 2 at the centres of cool-core clu sters to ~10 4 in the ou termost 
regions of the ICM (for a review, see lCarilli & Ta ylor 2002). Force 
balance then implies 

dlnp g 



dz 



(10) 



so that the inverse of the ion sound-crossing time across a thermal- 
pressure scale-height is equal to the dynamical frequency: 



Wdyn = [ 



1/2 



Vth 

H ■ 



(11) 



In general 6-VT 7^ 0, and so there may be a heat flux in 
the background state. In order to ensure our background state is in 

2 The numerical prefactor in equation (|5} depends on the exact form 
of the collision operator used. A more precise numerical prefactor is 
3075/1068 ~ 2.88 (e.g. ICatto & Simakovlliooll) : the ion collision fre- 
quency is Vi = 4y / 7rrije 4 Ai/3m i 1 '' 2 (/cbT;) 3 / 2 , where A; is the ion 
Coulomb logarithm. 

3 This assumes equal ion and electron temperatures, an assumption that 
may not hold in the outermost regions of galaxy clusters. 
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equilibrium, we must formally assume b-VQ = 00 However, as 
long as the timescale for the evolution of the background (global) 
state is longer than the local dynamical time, our results do not 
depend critically on the system actually being in global steady state. 

We allow perturbations (denoted by a S) about the background 
state and order their amplitudes as follows: 

Sv Sp ST 1 8p 1 SB 



■M, 



(12) 



v th p i M p /3 1 / 2 b /3V2 

where M <C 1 is the Mach number. This amounts to the Boussi- 
nesq approximation (i.e. relative changes in the pressure are much 
smaller than relative changes in the temperature or density). Sound 
waves are then eliminated from the analysis and so the flow be- 
haves as though it were incompressible, a good approximation in 
the ICM where typical velocities are much smaller than the sound 
speed. 

The perturbations are taken to have space-time dependence 
exp(<ri + ifc-r), where the growth rate a may be complex and the 
wavevector k — k x x + k y y + k z z. We order the timescales and 
the spatial scales as follows: 



The linearised entropy equation d 1 3t deserves special atten- 
tion. The first term on the right-hand side is responsible for the 
MTI. If the temperature increases in the direction of gravity, any 
alignment between the perturbed magnetic field direction and the 
temperature gradient (8b z 7^ 0) is unstable as long as conduction is 
rapid enough to ensure approximately isothermal field lines. When 
the temperature decreases in the direction of gravity, this term is 
stabilising. The second term on the right-hand side is responsi- 
ble for the HBI. If the temperature decreases in the direction of 
gravity, convergence/divergence of heat-flux-channeling magnetic 
field lines (k-Sb 7^ 0) leads to buoyantly-unstable density pertur- 
bations. When the temperature increases in the direction of gravity, 
this term is stabilising . 

Using equations l |20t and il9\ , the linearised internal energy 
equation d 1 8b becomes 



(<T + CJcond 



Sp 



Sv z — 4— In 



5 dz 



5/3 



SB Z - 2b z SB {l 1 dlnT 



B 



dz 



(21) 



Wdyn ~ kJcond 
1 M 



MH 



\mfp 



^visc ~ kvA ~ M kvth ~ M vu (13) 

(14) 



The latter ordering means that the relevant wavelengths are inter- 
mediate between micro- and macroscopic, viz. fc(A m f p Jf) 1 ' 2 ~ 1. 
Note that we are formally treating (m c /mi) 1//2 as a parameter of 
order unity here - a subsidiary expansion with respect to it will be 
done later. This completes the formulation of the problem. 



3.3 Linearised equations 

With account taken of the ordering introduced in Section 13.21 the 
linearised versions of equations ([T}h|3} and (|5}-(|7} are then 



k-Sv = 0, 



tSv 



. 2 fS P ± 1 SB H \ 2 8B 



. 8 P 
P 



f 3fcf|« t 2 h 



Sv, 



aSB — ik^BSv, 



(15) 



(16) 
(17) 



a 5v z - —In ^-jtt = — k- [bSQ - Sbxcb z — ) , (18) 



5 dz 



5/3 



5p 



dz 



fir 



SQ = -Xe 5b z Xc ifci 1 ST, 

dz 



Sp 
P 



5T_ 



(19) 



(20) 



where the subscript 1 1 denotes the vector component parallel to the 
background magnetic field (e.g. k\\ — b-k) and 8b = 8B±/B 
is the perturbation of the unit vector b. Equation d20| > expresses 
pressure balance for the perturbations, a consequence of our low- 
Mach-number ordering (eq. 112) . The total perpendicular pressure 
perturbation Sp± in equation l |16t is found by enforcing incom- 
pressibility (eq.1151. 



to leading order in M. Here we have introduced the characteristic 
conduction frequency, 



2 , 2 X cT 3.2 As 



^cond = — - Jq j- ( g-^- ) fc||A mfp _ffo; dyn 



nil 



1/2 



10 fc|| Amfpi? Wdyn, 



(22) 



where Ai (A c ) is the Coulomb logarithm of t he ions (electron s). 

Equations {[5)-l|2lJ differ from those in Ouataertl d2008l) only 
by the final term in the momentum equation d 1 6b . which is due 
to the perturbed pressure anisotropy (Braginskii viscosity). This 
term introduces a characteristic frequency associated with viscous 
damping: 

_ 3 fci|fth 3 2 

^visc = — — — = — Amfp/f Wdyn, (23) 

which is a factor ~6 smaller than uj cou d- 



3.4 Dispersion Relation 

The dispersion relation that results after combining equations dl5t - 
d 1 7b and d 3 1 b may be written in the following form: 



k 2 



a 2 (a + u cond ) + aN 2 



k 2 



dlnT K,' 

dz k 2 



2 b'xky 



a (a + u) cond ) + aN hWcondJ 



dlnTb 2 x k 2 



dz k 2 



(24) 



where k\ = k 2 — fc 2 is the square of the wavevector component 
perpendicular to the background magnetic fieldjfl 



~2 _ 2 , , 2 2 

a = a +K||«a ; 



(25) 



4 Another approach (see Balbus & Reynolds 2010) is to construct an equi- 
librium state in which conductive heating is balanced by radiative cooling. 



5 In c ontrast to our notation, Ouataert d2008l) and iBalbus & Reynolds] 
use k 2 ^ to denote the square of the wavevector component perpen- 
dicular to gravity, not to the background magnetic field. 
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and 



K = (1 -2b 2 z ) (k 2 x + k' 2 y ) + 2b x b z k x k z 



b x k — fcx + b x k y 



-b 2 z k 2 + kfi + biky. (26) 



We have written K, in three equivalent forms, all of which will prove 
useful in our analysis. We have also introduced the Brunt-Vaisala 
frequency given by 



(27) 



Were conduction, pressure anisotropy, and the magnetic field all 
to be ignored, equation d24t would reduce to the usual dispersion 
relation for internal gravity waves, a 2 — ~N 2 (k' x + k'y)/k 2 . 

It will also be beneficial to have the equations for the pertur- 
bations at hand, written in the limit of fast conduction (oj con d S> 
Wdyn ~ c): 

Sp _ ST ^ dlnT 
7" ~ ~~ ~ ^ dz 



a 2 (k z — 2fe z fc||) - (Taj v i sc b x (bxk) y 
a 2 k z — (ja; v i S c b x (bxk) y + g 2b x b z k x 



(28) 



dlnT, 



a k x + (TCJvisc b z (bxk) y + g— t— (1 — 2b z )k. 

a 2 k Z ~ (JWvisc b X (bXk)y + g ^ 26a; b Z fcj; 



(29) 



— =ife||^~ife||^ — 



cr 2 fcy + advise k 2 + g ^ (1 - 2b 2 )k 2 



a 2 k z — crwvisc b x (bxk) y + g — - — 26^ b z fcz 
-2,2 ,2 dlnT 

(7 fc + (TOJvisc Kl + <? ; /C 



dz 



cr 2 fc z — advise b x (bxk) y + g — ■ 2b x b z k x 



dz 



SB Z 
B 



(30) 



(31) 



— = ifc||Cn ~ ifc N 6 



a k\\+g 



dlnT 

dz 



bx k x 



a 2 k z — (ja; v i S c b x (bxk) y + g ^f"^ 26^ 6 Z fc^ 



dz 



(32) 



— = ifc N £ ± ~-^Ly + (yxb) 



cf (bxk) y + <jLj v i BC (bxk) y — g 



dlnT 

dz 



6 - A' r 



cr 2 fc z — advise b x (bxk) y + g—r—2b x b z .k x 



(33) 



where £ = cr<5u is the Lagrangian displacement of a fluid element 
(£ z > is upward). Equations d28t and ((32} imply that the La- 



grangian change in the temperature of a fluid element is 

dlnT _ 2i6 z dlnT SB ]{ 



AT ST , dlnT m 



dz 



•I 



dz B 



(34) 



emphasizing that perturbations in magnetic-field strength go hand- 
in-hand with changes in temperature. This will turn out to be an 
extremely important property for understanding the results of Sec- 
tion!] 



3.5 Nature of perturbations 

If we set ctMsc = 0, equation J24b re turns the standard MTI-HBI 
dispersion relation (see eq. 13 of louataertll2008b : 



~2 / s . AT 2 k x + k y 
a (a + Wcond ) + cr N — 



+ Wcond g- 



dlnT K 



dz k 2 



0. 



(35) 



The a = branch of this dispersion relation represents Alfven 
waves that are polarised with SB along the y-axis. They are un- 
affected by buoyancy. The other three modes are coupled slow 
and entropy modes. If we further take the limit of fast conduction 
(w<xmd w dyn ~ cr), the entropy mode becomes a ~ 
while the slow modes satisfy 



dlnT K 

dz k 2 



(36) 



When the temperature gradient and /C have opposite signs, one of 
the slow modes may become unstable to the MTI/HBI. 

Braginskii viscosity modifies this picture in two ways. First 
consider the limit k y — 0, in which the wavevector lies entirely in 
the plane spanned by gravity and the background magnetic field. In 
this case, equation d24t becomes 



a 2 (cr + Wcond) + cri^visc-^- (cr + W C ond) 



+ aN 



2 H~ fc'y 

~~k 2 



+ <^cond g- 



dlnT K 



dz k 2 



= 0. 



(37) 



The Alfven-wave branch of the dispersion relation is unchanged, 
since Braginskii viscosity does not affect motions perpendicular to 
the magnetic field. In contrast, slow-mode-polarised perturbations 
with kj_ ^ are damped. This property of Braginskii viscosity is 
the root cause of the significant changes to the nature of the allowed 
unstable HBI (Section l4~TTTT > and MTI (Section l42~T1 l modes. 

Next consider the general dispersion relation (124) with k y 7^ 
0. In this case, Braginskii viscosity couples the Alfven and slow 
modes. This can be seen most clearly by taking the fast-conduction 
limit (oj con d 2> tJdyn ~ cr) of equation J24t : 



-2I~ 2 k 2 ± dlnT K. 



dlnTb 2 x k 2 



dz k 2 



(38) 



When cj v i sc ^> ojd yn the slow mode is rapidly damped, leaving 
only crojvisc k\jk 2 to highest order in the parentheses on the left- 
hand side of equation d38b . This term cancels the similar factor on 
the right-hand side, ultimately leading to 



dlnT b 2 x k 2 y 
'~d~z W' 



(39) 
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which may be unstable when the temperature increases in the di- 
rection of gravityQ We will elaborate on this result in Section l4.2.2l 
where we discuss this new Alfvenic' version of the MTI, but for 
now we explain the physical content of equations d38fc and l |39l >. By 
damping motions along field lines, Braginskii viscosity effectively 
reorients magnetic field perturbations to be nearly perpendicular 
to the background magnetic field (via flux freezing). These modes 
therefore display characteristics of both slow and Alfven modes: 
they have density and temperature perturbations, and therefore are 
subject to buoyancy forces, but their velocity and magnetic field 
perturbations are predominantly polarised across the mean field. 

We note in passing the striking similarity between equation 
l |38t and the dispersion relation for the axisymmetr ic magnetoro- 
tational instability su bject to Braginskii stresses (Balbusl liooi ; 
llslam & Balbusll2005h ; 

.2 fci dlnfi 2 k% 



HBI without Braginskii 



+ advise 



A- 



+ g- 



dR k 2 



dm ft 2 bik 



:9- 



dR k 2 



(40) 



where g — Q 2 R in a rotating disc. Aside from a AO. 2 tenn due 
to epicyclic motions, the equivalence is revealed by relabelling 
the disc coordinate system (R, (f>, Z) -f-> (z, x, y) and swapping 
one free energy source (temperature gradient) for another (angu- 
lar velocity gradient) Q Braginskii viscosity couples the Alfven- 
and slow-mode branches of the dispersion relation via the an- 
gular velocity gradient in very much the same way that it cou- 
pled these branches via the temperature gradient in equation | |38I >. 
Furthermore, when the angular velocity decreases outwards and 
i^visc S> ^ ~ <J, the first term on the right-hand side of equa- 
tion d40b due to Braginskii viscosity overwhelms the second term 
due to epicyclic coupling and drives the magnetoviscous instability 
(MVI) by endowing slow-mode perturbations with Alfven-mode 
characteristics (as in eq. 138) . We therefore identify the behaviour 
revealed by equation i39\ as the temperature-gradient analog of the 
MVI. 



4 RESULTS 

4.1 dT/dz > 0: Heat-flux-driven Buoyancy Instability 

We first investigate the effects of Braginskii viscosity on the stabil- 
ity of a stratified atmosphere in which the temperature decreases in 
the direct ion of gravity, i. e. dT/dz > 0. Such an atmosphere was 
shown bv lOuataeril d2008l) to be susceptible to the HBI if K. < 0. 

4.1.1 Case ofk y = 0: Standard HBI with and without Braginskii 
viscosity 

If Braginskii viscosity is ignored, it is straightforward to show from 
equation ([35) that the maximum HBI growth rate 



6 There are other instances of an anisotropic damping mechanism coupling 
the Alfven and slow mode branches of a dispersion relation via a free energy 
gradient. In weakly-ionised plasmas, the interaction between velocity shear 
and anisotropic magnetic resistivity (ambipolar diffusion and the Hall ef- 
fect) results in such a cou pling - one which ultimately leads to shear-driven 
instabilities <Kunzl200Sl) 

7 We refer the reader to balbusl fcOOfj 1200 lb for a cogent discussion of 
the analogy between angular momentum and entropy that underlies these 
mathematical similarities. 
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HBI with Braginskii 
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Figure 3. HBI growth rate (normalised to the maximum growth rate 
^/ g d In T/dz) without (top) and with (bottom) Braginskii viscosity for 
a stratified thermal layer with dlnT/dlnp = — 1 threaded by a vertical 
magnetic field (b z = 1). Magnetic tension is neglected; its effect is dis- 
cussed in Section 14.1.31 Each contour represents an increase in the growth 
rate by 5 per cent. The solid lines (given asymptotically by eq. |52| with 
Braginskii viscosity and eq. l43l without Braginskii viscosity) trace the max- 
imum growth rate for a given total wavenumber k; the maximum growth 
rate is given by eq. 1491 with Braginskii viscosity and asymptotically at 
k|| (^mipH) 1 ^ 2 3> 1 by eq. 1411 without Braginskii viscosity. Braginskii 
viscosity dramatically reduces growth rates everywhere except for a narrow 
band of wavenumbers around fci | given by eq.f52] Galaxy clusters with cool 
cores typically have H/A m f p ~ 10 2 -10 3 at radii for which the tempera- 
ture profile increases outwards, and so the maximum wavenumbers for each 
axis span the range kH ~ 40-127. 



°"HBI.: 



dlnT, 2 



occurs for wavevectors satisfying 

fcf, 



, 2 CHBI.max 
9z 

^cond 



1 + 



1 d lnp 



5 dlnT 



< 1 



(41) 



(42) 



k\ ~ 

to leading order in Wdyn/ucond, where we have assumed k\\H <g 
bz/3 1//2 - i.e. magnetic tension is negligible on the scales of inter- 
est. Equation d42b reveals that the HBI has a strong preference for 
perpendicular wavenumbers. More precisely, using the definition of 
u coll( i (eq.l22t in equation J42t . we find that the maximum growth 
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rate occurs along a path through fc-space on which 



x 0.6fc? /4 ( 1 + 



dlnp 



dlnT 



1/4 



dlnT 



1/8 



d hip 



(43) 



This behaviour is exhibited in the top panel of Fig. [3] which shows 
HBI growth rates in the (fen, k±) plane for b z = 1 (without Bra- 
ginskii viscosity). The solid line in the plot traces the maximum 
growth rate through wavenumber space; it quickly asymptotes to 
equation d43 1 . 

The HBFs preference for perpendicular wavenumbers is also 
reflected in the corresponding eigenvectors. Using equations d28t 
and d36t we find that the density perturbation associated with the 
HBI, 



5p 
P 



dlnT k 2 

dz k'i 



M 

A- 2 



(44) 



is greatest when fc 2 <C k 2 so that, e.g., upwardly-displaced fluid 
elements have the largest possible decrease in their density. More- 
over, perpendicular wavenumbers are necessary to generate linear 
perturbations in magnetic field strength, 

8Bu 



3 I 



B 



k' x 



(45) 



which lead to local convergence/divergence of the background heat 
flux and consequent heating/cooling of the plasma (eq.|34|). 

The problem is that it is precisely such perturbations that are 
damped by Braginskii viscosity (see the bottom panel of Fig. [5}. By 
equation d34t , upward displacements along magnetic field lines go 
hand-in-hand with local heating (AT > 0) and a local increase in 
the magnetic field strength (8B\ | > 0)0 This causes a negative vis- 
cous stress that damps motions along field lines, thereby rarefying 
the magnetic field and reducing the strength of the perturbed heat 
flux. This can be seen quantitatively by explicitly writing down the 
buoyancy and viscous forces in the z-component of the momentum 
equation d!6t for the simple case b z — 1: 



d 2 ^ _ , dlnT. d£ z 



(46) 



For wavenumbers satisfying the ordering cj v i S c 3> i^dyn, or 
fc^Amfp-ff S> 1, it is straightforward to show from equation (1381 
that the growth rate is 



g dlnT 

Jvisc dz 



2 

^dyn 



(47) 



to leading order in ojdyn/Wiso In other words, the buoyancy and 
viscous forces become nearly equal and opposite as the plasma be- 
comes more and more collisionless. The growth rate decreases ac- 
cordingly. 

Despite all this, there are modes that remain unstable to the 
HBI and retain non-negligible growth rates. However, it turns out 
that they are confined to a thin band of wavenumber space in which 
conduction is fast but viscous damping is small: 



^cond ^ ^dyn ^ ^visc , 



or, using the definitions d22t and d23t . 

3fc||(A mfp #) 1/2 > 1 > fc||(A mfp #) 1/2 . 



(48a) 



(48b) 



8 While Alfvenically-polarised modes suffer no viscous damping, they are 
HBI stable because 6B\ i = 0. 



Using the fact that k± ~ k for the fastest-growing Braginskii- 
HBI modes, it is possible to obtain analytic solutions for the max- 
imum growth rate and fastest growing wavenumber. Defining e = 
ixMsc/Wond ~ 0.1, the maximum growth rate 



C"HBI,max 



1 - 2e 



1/2 



2 dlnp 
5 dlnT 



2e d In p 

!1+£ + y diiTr 



1/2 



1/2 



d lnp 



dlnT 



occurs at a parallel wavenumber satisfying 

c(l + e) 



£| | ^mf'p 



H 



2 2 
°"HBI,max — °mi 



-1/2 



3 0" ma x UJdyn 



dlnT 



d lnp 



1/2 



O(e) 



1/2 

(49) 

(50) 
(51) 



Since k\\ — k x b x + k z b z , equation \5i\ implies that the maximum 
growth rate is attained along two straight lines in the (k x ,k z ) plane 
given by 



fc z (A m f p -f/') 1/ ' 2 ~ t~ k x {\ mlp H) 1/2 



±£ 



1/4 



1/2 



dlnT 



dlnp 



1/4 



+ 0(e 



3/4x 



(52) 



This behaviour can be seen in Fig. [4] which exhibits HBI growth 
rates in the (k x ,k z ) plane with (upper row) and without (lower 
row) Braginskii viscosity for k y — and various magnetic field ori- 
entations. The solid lines trace the maximum growth rate through 
wavenumber space; they quickly asymptote to equation d43t with- 
out Braginskii viscosity and equation d52t with Braginskii viscos- 
ity. 

For a fiducial cool-core temperature profile dlnT/dlnp = 
— 1, equations d49t and d50l > give a max = 0.576 2 uJdyn and 
k\\ (A m fp7T) 1//2 = 0.606?/ 2 , respectively. With typical values 
of H/Xmip ~ 10 2 -10 3 in the inner ~200 kpc of cool-core 
clusters where the temperature increases with height, this im- 
plies fc||7T ~ 6b z //2 -196 z //2 (increasing inwards). These modes 
are quite extended along the magnetic field direction and can- 
not be considered local. This is likely to have important impli- 
cations for the non-linear evolution of the HBI, particularly as 
the HBI reorients the mean magnetic field to be more and more 
perpendicular to the temperature gradient. For example, taking 
dlnT/dlnp = —1 and /T/A m f p = 200, the parallel wavelength 
of maximum growth A|| max is equal to the thermal-pressure scale- 
height H when the magnetic field makes an angle of 6 ~ 33° 
with respect to the x-axis. Thus, the field-line insulation found 
by many num erical simulations to be a consequence of the stan- 
dard H BI (e.g. lParrish. Ouataert & Sharmal2009tlBogdanovic et all 
2009) may not be as complete as is currently believed. 

Note further that equation ( 151) in the limit e — > does not re- 
duce to the no-Braginskii case (eq.|43|l. Moreover, the relationship 
between fc| | and k± for the fastest growing modes discontinuously 
changes from k^ oc k 1 / 2 without Braginskii viscosity (eg. 143 \ to 



k 



const with Braginskii viscosity (eq. 151b for perpendicular 



wavenumbers satisfying 



k± (A mfp /T) 1/2 > 0.5 



dlnT 



i/i 



dlnp 



a/2' 



(53) 
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Figure 4. HBI growth rate (normalised to yj g d InT/dz) for fc^ = 0, dlnT/dlnp = —1, and various magnetic field orientations 6 = cos^ 1 (b x ). 
Magnetic tension is neglected; its effect is discussed in Section 14.1.31 Braginskii viscosity is included in the bottom row of plots. Each contour represents an 
increase in the growth rate by 5 per cent. The dashed line denotes the direction of the background magnetic field. The solid lines (given asymptotically by eq. 
[52] with Braginskii viscosity and eq. [43] without Braginskii viscosity) trace the maximum growth rate for a given total wavenumber k; the maximum growth 
rate is given by eq.|49]with Braginskii viscosity and asymptotically at fc|| (^mfpH) 1 ^ 2 2> 1 by eg. 14 II without Braginskii viscosity. The overall maximum 
growth rates in each of the Braginskii-HBI plots (bottom row) are reduced by a factor ~1.66 relative to those in the respective standard HBI plots (top row). 
The thick black diagonal region in each plot is where fen =0. 



This reflects the fact that including fast anisotropic heat conduc- 
tion while neglecting Braginskii viscosity is a singular limit of the 
equations. 



and a pure oscillation if the left- and right-hand sides of these in- 
equalities are in fact equal. Indeed, in the limit Wdyn/w-visc <C 1 
equations J28 M341 imply 



4.1.2 Case ofb x k y / 0: Alfvenic HBI 
If b x k y 7^ 0, the situation is actually worse: 

dlnT 



dlnT b 2 x k 2 y \ 1/2 u 2 d 
°" — 1 I 9— I + 



dz k 2 



d lnp 



2k\ 



Ir k 

U x Hjy 



K_ 

k 2 
(54) 



to leading order in cjdyn/Wisc <C 1. The HBI becomes a slowly- 
growing overstability for wavevectors satisfying 



ky > 



bx 



fcii(fcX6)j 



(kXbfy 



a weakly-damped oscillation for wavevectors satisfying 



k y < 



-(kxbf y 



(55) 



(56) 



P dz \Wvisc 



+ o ' Udyn 



AT 
~~f~ 

5B \\ 
~B~ 



iO 



+ o 



2 

w dyn 
^visc 



(57) 



(58) 



(59) 



By effectively reorienting magnetic field perturbations to be nearly 
perpendicular to the background magnetic field, Braginskii vis- 
cosity prevents slow-mode perturbations from tapping into the 
free energy carried by the background heat flux. To highest or- 
der 111 LOdyn /ajvisc, these modes appear as HBI-stable Alfven waves 
whose magnetic tension has been effectively increased by the ad- 
verse temperature gradient. 
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4.1.3 Effect of magnetic tension on the HBI 

When both Braginskii viscosity and magnetic tension are included 
there are two parallel-wavenumber cutoffs, the relative magnitude 
of which may play an important role in the evolution and non- 
linear saturation of the HBI. Roughly speaking, magnetic tension 
is appreciable for parallel wavenumbers satisfying fc|| imax ^A > 



MTI without Braginskii 



o"max, where fen max and <j ma x are given by equations d50t and 
d49t respectively. For a fiducial cool-core temperature profile, 
dlnT/dlnp = —1, this amounts to an upper limit on the plasma 
beta parameter of b z j3 < H/X m f p . For /3 less than this value, mag- 
netic tension - not Braginskii viscosity - sets the fastest-growing 
mode. 



4.2 dT/dz < 0: Magnetothermal Instability 

Next we investigate the effects of Braginskii viscosity on the stabil- 
ity of a stratified atmosphere in which the temperature increases in 
the direct ion of g r avity, i.e. dT/dz < 0. Such an atmosphere was 
shown by BalbuJ ( feOOOh to be susceptible to the MTI if K, > 0. 



4.2.1 Case of k y = 0.' Standard MTI with and without 
Braginskii viscosity 

Consider first the case k y = 0. Equation J3 8 b shows that the maxi- 



mum MTI growth rate 



. Ilia 



9 



dlnT 



dz 



(60) 



,1/2 



occurs for k± = 0, where we have assumed that k^H <C b x /3 
(i.e. magnetic tension is negligible on the scales of interest). The 
physical reasons for this are simple. Since k-8B — 0, taking k± — 
implies SB^ = 0. This not only ensures that any background 
heat flux is unable to cool (heat) upwardly (downwardly) displaced 
fluid elements (see eq,l34t. but also that pressure anisotropy cannot 
damp these modes. 

If we allow for a small wavenumber component perpendicular 
to the background field, k\ <C k 2 , the leading-order solution for 
the growth rate is given by 



°"MTI. 



1 - 



A: 



b 2 k 2 



- CMTI. 



max ^visc 



TT- (61) 



The first negative contribution on the right-hand side of this equa- 
tion is tied to the fact that generating a <5B|| implies k-Sb ^ 0. 
Having k-5b ^ causes upwardly displaced fluid elements, which 
are trying to heat and rise by the MTI, to be cooled by the locally 
divergent background heat flux (and vice versa for downward dis- 
placements). This reduces the efficiency of the MTI and implies 
that unstable modes with large growth rates are confined to a wedge 
in wavenumber space of width 



k± < b x k 



II- 



(62) 



As fen increases, more and more perpendicular wavenumber space 
becomes available for fast-growing modes (see top panel of Fig.O. 

When Uvisc S> c^dyn, however, this term is relatively unim- 
portant when compared to the last term on the right-hand side of 
equation J6U . Modes with k\ 7^ are rapidly damped by Bragin- 
skii viscosity. This behaviour continues beyond the small values of 
k± , where the growth rate then becomes 



9 



dlnT (btk 



dz 



A' 



- 1 



(63) 
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Figure 5. MTI growth rate (normalised to the maximum growth rate 
\f—Q d lnT/dz) for k y = without (top) and with (bottom) Braginskii 
viscosity for a stratified thermal layer with dlnT/dlnp = 1/3 threaded 
by a horizontal magnetic field (b x = 1). Magnetic tension is neglected; its 
effect is discussed in Section l4.2.3l Each contour represents an increase in 
the growth rate by 5 per cent. Braginskii viscosity suppresses MTI growth 
rates for k 7^ fen (see eqs 1611 and 1641 . Galaxy clusters typically have 
H/X m fp ~ 10-100 at radii for which the temperature profile decreases 
outwards, and so the maximum wavenumbers for each axis span the range 
kH ~ 13-40. 



to leading order in ajdyn/^visc 1. These effects are evident in 
Figs [5] and [6] which exhibit growth rates in the (fc^, k z ) plane for 
various inclinations of the background magnetic field. From equa- 
tion l |63t , we find that unstable modes with large growth rates are 
confined by Braginskii viscosity to a narrow band in wavenumber 
space of width 



fex (A mfp H) 1/2 



dlnT 
dlnp 



1/4 



bV 2 



(64) 



about the line k\\b. In contrast with the no-Braginskii case (eq. 
|62t , going to larger fen does not open up more perpendicular 
wavenumber space for fast-growing modes. 
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Figure 6. MTI growth rate (normalised to y/ — g dlnT/dz) for k y = 0, dlnT/dlnp = 1/3, and various magnetic field orientations 9 = cos -1 (b x ). 
Magnetic tension is neglected; its effect is discussed in Section l4.2.3l Braginskii viscosity is included in the bottom row of plots. Each contour represents an 
increase in the growth rate by 5 per cent. The dashed line denotes the direction of the background magnetic field, which also traces the maximum growth rate 
for a given total wavenumber k. The maximum growth rate is given by equation J60t and occurs for kj_ = 0. Braginskii viscosity suppresses MTI growth 
rates for k ^ ku (see eqs l61l and l64t . 



4.2.2 Case ofb x k y / 0: Alfvenic MTI 
Recall equation d38b : 

dlnT K 



a + acj v i sc — + g- 



-auj v i sc g- 



dlnT b 2 x k 2 



dz 



(65) 



which is the general dispersion relation 1 1241 written in the fast- 
conduction limit (ojcond ^> Wdyn ~ c). When b 2 k 2 ^ 0, the right- 
hand side of this equation becomes active and leads to behaviour 
otherwise absent without Braginskii viscosity. The Alfven-mode 
branch of the dispersion relation is now coupled to the slow-mode 
branch; slow-mode perturbations induce an Alfvenic response. 
Consider further the limit o; v isc 2> Wdyn- Then we obtain eq. 

09): 



dlnT 



dz 



b 2 k 2 



(66) 



This is always unstable, regardless of the sign ofK,, and is maximal 
when ky = k\ (i.e. the projection of k onto the x-z plane is par- 
allel to the magnetic field). In other words, wavevectors with large 
parallel components and any component along the y-axis grow at 
CTMTi.max (see Fig.[7]l. Indeed, one can readily show from equation 
d65t that the growth rate for these ky = k\ Alfvenic MTI modes is 



CMTLmax 



dlnT fa 2 



a; v i sc dlnp 2 
to leading order in Wdyn/uVisc. so that one only requires 

'dlnT\ 1/4 b z 



fc||(A mfpJ fT) 1/2 > 



dlnp 



& l/2 



(67) 



(68) 



to bring the growth rate close to CMTi.max (e.g. see the rightmost 
panel in the bottom row of Fig. 0. One consequence is that it is 
no longer necessary to go to very large parallel wavenumbers (and 
risk stabilisation by magnetic tension) just to marginally destabilise 
k y ^ modes at small b x (e.g. see the rightmost panel in the top 
row of Fig. [7}. 

The physical origin of this new behaviour may be uncovered 
by computing the eigenvectors d28t . &32\ and d34t to leading order 

in ojdyn /^visc : 



p 



dlnT 



5Bh 
B 



dz 

I^dyn 



1 - 



^dyn 



dlnT 
d hip 



dlnT 
d hip 



1/2 



AT ^ 

rp — _ V 



dlnT w d yn / dlnT 
dz td visc \ dlnp 



1/2 



1/2 2b| 
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(69) 
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Figure 7. MTI growth rate (normalised to yj— g A hyT/dz) for k y = k 2 ^, dlnT/dlnp = 1/3, and various magnetic field orientations 9 = cos -1 (b x ). 
Magnetic tension is neglected; Braginskii viscosity is included in the bottom row of plots. Each contour represents an increase in the growth rate by 5 per cent. 
Many k y ^ modes that are either stable or only grow slowly in the absence of Braginskii viscosity become unstable with growth rate <tmti max when 
Braginskii viscosity is included. 



These are essentially MTI modes that have been freed from the 
unfavourable consequences of having k-Sb ^ by rapid paral- 
lel viscous damping. Rapid Braginskii damping endows buoyantly- 
unstable slow-mode perturbations (eq.|69t with Alfvenic character- 
istics (i.e. perturbed magnetic fields and velocities that are predom- 
inantly oriented perpendicular to the background magnetic field - 
eq. 170) , which allows fluid elements to approximately maintain 
their temperature as they are displaced (eq. 171) . The presence of 
a non-zero k y is necessary in order to ensure 5B\ | is vanishingly 
small for arbitrary k x and k z , while simultaneously preserving the 
divergence-free constraint on the magnetic field. Therefore, many 
wavevectors for which K, < that are stable to the standard MTI 
(e.g. ky — k\ > k 2 b x /b 2 z , which is the black region in the 
upper-right panel of Fig. [7} are actually unstable with growth rates 

— C"MTI,max- 



4.2.3 Effect of magnetic tension on the MTI 

When k y — the fastest growing MTI modes, those with k = k\ i , 
are unaffected by Braginskii viscosity. In this case, magnetic ten- 
sion provides the only parallel-wavenumber cutoff by suppressing 



wavenumbers for which k\ 



VA 



kH = kuH > p 



1/2 



dlnT 
d hip 



CTMTI,] 



1/2 



(72) 



Setting kH ~ 2-k provides a strict lower limit on beta below which 



dlnp / 2n 
dlnT I 67 



(73) 



For a fiducial cluster temperature profile beyond the cooling radius, 
dlnT/dlnp = 1/3, this gives f3 < 120 b x 2 . 

If we relax our restriction on k y , Braginskii viscosity drives 
Alfvenic modes unstable by coupling them to buoyantly-unstable, 
rapidly-damped slow modes. Setting knVA ~ fMTl.max in equa- 
tion J68t we find that, unless 



H 



d lnp 



1/2 



(74) 



A mfp \dlnT, 

there are fast-growing Alfvenic MTI modes. With typical values of 
H/A mf p ~ 10-100 and /3 ~ 10 3 -10 4 in the outer regions of the 
ICM where the temperature decreases outwards, magnetic tension 
is unlikely to affect these modes except possibly when the field is 
nearly vertical (b x 1). Note that increasing k y does not increase 
magnetic tension since b _L y. 

4.3 dT/dz < 0: Heat-flux-driven Buoyancy Overstability 

The reader may have noticed a very thin vertical band of unstable 
modes in the 9 — 60° panels of Fig. [7] These modes, found recently 
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bv lBalbus & Reynolds! fioiol) . are g-modes driven overstable by a 
background heat flux. They become important only when the back- 
ground magnetic field is vertical {b z = 1, K = — kjj, since this 
field orientation is (linearly) stable to the MTI. In this Section we 
analyse these modes while including Braginskii viscositylf] 

We begin by finding the fastest-growing mode in the absence 
of Braginskii viscosity. Our task is greatly simplified by knowing 
a priori that k± ~ k for these modes. We also neglect magnetic 
tension; we will verify a posteriori that it is unlikely to affect the 
fastest growing mode for conditions found in the outer regions of 
the ICM where this overstability may be present. Our dispersion 
relation d24t then becomes 



<t 3 + a 2 uj con d + aN 2 



■ <^cond g 



dlnT 



dz 



0. 



(75) 



Solutions of this equation have both real and imaginary parts, 
a = 7 + iuj. Substituting this decomposition of a into equation 
l |75t and separating into real and imaginary parts gives two equa- 
tions for 7 and lo as functions of oj con d. The maximum growth rate 
is then found by differentiating these with respect to u; con d, setting 
d~y I 'Second = (so as to maximize the real part of a), and solv- 
ing these four equations simultaneously. Here we simply state the 
resultEB 

1/2 



7max ■ — k^dyn 



1 /dlnT 



d lnp 



dlnp 
dlnT 



ldlnT 1 dlnp 



10 4 dlnp 100 dlnT 
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dlnp 

^cond.max — ^dyn T" I "T"; ™ 

5 V d In T 
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so that 



''I | ,max 



(XmlpH) 



1/2 



0.1 



dlnp 
dlnT 



1/4 



(76) 



(77) 



(78) 



(79) 



(80) 



This mode is overstable if 
dlnT 1 
d lnp 5 

In order for magnetic tension to significantly affect this mode, 
k\ | ,max«A > 7max or, equivalently, 

P < 0.08 



H /dlnT 1 V 2 /dlnT^ 1/2 



d lnp 



(81) 



A m fp V dlnp 5, 

With typical values of #/A mfp ~ 10-100 and /3 ~ 10 3 -10 4 in 
the outer regions of the ICM where the temperature decreases out- 
wards, it is highly unlikely magnetic tension will affect the growth 
rate of this mode. 

If we take into account Braginskii viscosity, it is straightfor- 
ward to show that a necessary condition for stability is given by 

(82) 



r , / htasc kj_ htasc i ' „ i 



where 



9 We have chosen not to present a similar analysis for the radiative- 
cooling-driven g-mode overstability also found by Balb us & Revnoldsl 
12010) for dT/dz > 0. The necessary assumption that the local cooling 
rate is comparable to the dynamical frequency conflicts with our assumption 
of an equilibrium background state that evolves slower than the instabilities 
do. 

10 This assumes N 2 > or, equivalently, dlnT/dlnp < 2/5. 



,'2 dlnp 
- 1 5 dlnT 



(83) 



If the temperature profile satisfies equation (180), then 1Z' < 
and b uoyant modes are potentially overstable (Balbus & Reynolds 
2010). A comparison with equation (27) of iBalbus & Revnoldsl 
fcoich reveals that Braginskii viscosity modifies this by effectively 
increasing 7Z' by ~cj v i sc /aj con d ~ 0.1 so that equation d80t be- 
comes 

I _ _ ZliZ I ~ 0.22 



dlnT 1 
dlnp 5 



2 CJ C ond + ^vi 

and by further stabilising modes for which 

WviscW con d > —N 2 1Z', 

or, using the definitions d22t and J23 1 . 



fc||(A mfpJ f7) 1/2 >0.6^ 



dlnT 



lnp 



1/4 



(84) 



(85) 



(86) 



However, the maximum growth rate of the overstability changes 
very little when Braginskii viscosity is included: 



7 — 7max — U>dy 



d lnp 
20 \ dlnT 



1/2 



i + 1 dln P 
5dlnT 



1 

25 



d lnp 
dlnT 



(87) 



to leading order /ucond ~ 0.1. For a fiducial cluster 

temperature profile of dlnT/dlnp = 1/3, the Braginskii term 
amounts to a correction <16 per cent. Braginskii viscosity does 
not significantly affect the fastest-growing overstable mode. 



5 DISCUSSION 

The low degree of collisionality found in astrophysical plasmas 
such as the ICM causes heat and momentum transport to be- 
come anisotropic with respect to the magnetic field direction. This 
implies anisotropic heat flux and pressure. The former has been 
previously found to play a destabilising role in thermally strati- 
fied atmospheres, causing instabilities such as the MTI (when the 
tempe rature increases in the direction of gravity; Balbusl I2000L 
120011) and the HBI (when the temperature decreases in the direc- 
tion of gravity; [Q uatae rj ||2008l) . as well as <?-mode overstabilities 
jBalbus & Revnoldsll2010h . In this paper we have concentrated on 
the consequences anisotropic pressure has for the stability of the 
ICM. 

We have argued that one cannot consider the limit of fast con- 
duction along field lines while neglecting the Braginskii pressure 
anisotropy. Although there is a timescale disparity between the two 
effects - anisotropic heat conduction acts on a timescale a factor 
~10 shorter than does pressure anisotropy - both are generally 
much faster than (or at least as fast as) the dynamical timescale. 
Since the MTI and HBI occur with a growth rate comparable to 
the dynamical frequency, pressure anisotropy affects their dynam- 
ics significantly. 

In the case of the HBI, its propensity (or, more accurately, 
its need) to generate fluctuations along the background magnetic 
field suffers from the requirement for particles in a weakly colli- 
sional plasma to conserve their first and second adiabatic invari- 
ants. The HBI changes the field strength to linear order, which in- 
duces a pressure anisotropy, which manifests itself as Braginskii 
viscosity and kills off the motions that generated the change in 
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field strength in the first place. The only motions to entirely es- 
cape this constraining effect of pressure anisotropy - those that are 
Alfvenically polarised - are also those that are stable to the HBI. 
The fastest growing HBI modes no longer occur at large parallel 
wavenumbers, but rather at wavenumbers satisfying the timescale 
ordering cj cond > u dyn > oj visc , or 

3fc||(A mfp tf) 1/2 > 1 > fc||(A mfpJ ff) 1/2 . 

Perturbations whose wavelengths along the background magnetic 
field are smaller than the thermal-pressure scale-height by at least 
a factor ~27r(A m fp/i/) 1//2 , while potentially unstable to the HBI, 
are nevertheless strongly damped. Small-wavelength perturbations 
whose wavevectors have a component perpendicular to both grav- 
ity and the background magnetic field behave like modified Alfven 
waves that are only slowly growing or decaying (depending on their 
exact wavevector orientation; see eq.|54). Unless b z ft < i//A m f p , 
Braginskii viscosity - not magnetic tension - sets the maximum 
unstable parallel wavenumber. 

The situation with the MTI is more complicated. The standard 
MTI has a slight preference for wavevectors with projections in 
the x-z plane that are aligned with the background magnetic field 
(see the top row of Fig. [6}. This obviates heat exchange with any 
background heat flux, a stabilising effect when the temperature in- 
creases in the direction of gravity. Pressure anisotropy reinforces 
this preference, since perturbations whose projected wavevectors 
are not perfectly aligned with the background magnetic field are 
subject to strong viscous damping (see bottom row of Fig. |6). We 
have also found that many modes that were considered MTI-stable 
[e.g. ky = k\ > fc 2 b 2 /& 2 ] or slowly growing become unstable 
in the presence of Braginskii viscosity and grow at the maximum 
possible rate for a given background magnetic field orientation (see 
Fig- EJ- This is because, when k-(bxg) 7^ 0, Braginskii viscosity 
couples the Alfven- and slow-mode branches of the dispersion re- 
lation, so that slow-mode perturbations excite a buoyantly-unstable 
Alfvenic response. By damping perturbations along magnetic field 
lines, pressure anisotropy frees these modes from the unfavourable 
consequences of having local field-line convergence/divergence. 

We anticipate that many of the results found by numerical sim- 
ulations of the MTI and HBI will change both quantitatively and 
qualitatively when the equations including both anisotropic heat 
and momentum transfer are implemented. This will likely have im- 
portant consequences for our understanding of the thermodynamic 
stability of the ICM. 

Depending on the degree of collisionality in the cool cores of 
galaxy clusters, the field-line insulation found in many simulations 
to be a consequence of the non-l i near evolution of the HB I (e.g. 
IParrish. Ouataert & Sharmal 120091; iBogdanovic et al.ll2009l) might 
be attenuated. This is because the large wavenumbers required to 
keep the HBI in action as the magnetic field becomes more and 
more horizontal are strongly suppressed by the pressure anisotropy 
they generate. Moreover, the wavenumbers at which the HBI sur- 
vives largely unsuppressed have parallel components too small to 
rigorously be considered local, especially as the HBI reorients the 
mean field to be horizontal (see eq. |50t . For a fiducial cool-core 
temperature profile dlnT/dlnp = —1, the parallel wavelength 
of maximum growth is equal to the thermal-pressure scale-height 
when b z ~ 110 X m f p /H. It is therefore tempting to speculate that, 
in the absence of strong turbulent stirring by an external agent, there 
exists a link between the degree of collisionality in cool cores and 
the mean direction of the magnetic field. 

In the outer regions of non-isothermal clusters the non-linear 
evolution of the MTI may be more vigorous than previously 



thought, since many modes classified as stable or slow-growing 
are actually maximally unstable. Moreover, the fact that Braginskii 
viscosity couples damped t, / slow modes with MTI-unstable 
Alfven modes, a feature not present in the standard MTI, may pro- 
foundly affect the non-linear evolution of the magnetic field. On the 
other hand, the nonlinear excitation of the MTI out of its linearly- 
stable end state (b x — 0), which is triggered by buoyantly-neutral 
horizontal motions (see lMcCourt et aLll2010l) . is unlikely to be af- 
fected by Braginskii viscosity. These motions occur perpendicular 
to the magnetic field (i.e. k± = 0) and are therefore undamped by 
Braginskii viscosity. In either case, the spectrum of unstable modes 
will certainly be different, not only due to the presence of a paral- 
lel viscous cutoff but also because pressure anisotropy significantly 
modifies the dependence of growth rate on wavenumber. 

The h eat-flux-driven buoyancy o verstability elucidated ana- 
lytically bv lBalbus & Reynolds] J20icl) and numerically by T. Bog- 
danovic (private communication) is not significantly affected by 
pressure anisotropy. This is because the overstability occurs at suf- 
ficiently small wavenumbers such that conduction (and therefore 
viscosity) is not overwhelming (see eqs,|79|and!87t. Braginskii vis- 
cosity only shifts the stability boundary slightly (eq.|82>. 

The nonlinear evolution of the MTI and HBI in the presence of 
anisotropic viscosity is, at least in principle, amenable to numerical 
simulation. The Athena code affords one promising venue, as it is 
already set up for the inclusion of both anisotropic conduction and 
anisotropic viscosity (J. Stone, private communication). In prac- 
tise, however, the implementation of pressure anisotropy into a nu- 
merical code is rather nuanced. If the pressure anisotropy exceeds 
~l//3, very fast microscale instabilities (e.g. firehose, mirror) can 
be triggered, which will grow rapidly at the grid scale and wreak 
havoc upon a simulation if left unchecked. Exactly how such insta- 
bilities nonlinearly saturate remains very much an open question 
(e^g. see lSharma et al ]|2006l;ISchekochihin et alj|200 8; Rosi n et al.l 
|2010|) and, in lieu of performing a full kinetic calculation, important 
choices will need to be made by the simulator regarding anisotropy 
limiters. Despite these rather foreboding complications, properly 
simulating the ICM with equations that include both anisotropic 
heat and momentum transfer would be a major step forward in our 
understanding of the dynamical stability of the ICM. 
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